Sensors 2012, 12, 7975-7993; doi:10.3390/sl20607975 



OPEN ACCESS 



sensors 

ISSN 1424-8220 

www.mdpi.com/journal/sensors 

Article 

Reliable Finite Frequency Filter Design for Networked Control 
Systems with Sensor Faults 

He-Hua Ju \ Yue Long ^ and Heng Wang 

^ College of Electronic Information and Control Engineering, Beijing University of Technology, 

Beijing 100124, China; E-Mail: juhehua@bjut.edu.cn 
^ College of Information Science and Engineering, Northeastern University, Shenyang, Liaoning 

110004, China; E-Mail: lemonshell@yahoo.cn 

* Author to whom correspondence should be addressed; E-Mail: wangheng@bjut.edu.cn or 
wanghengl9801980@126.com; Tel.: -1-86-1581-1005-645. 

Received: 6 April 2012; in revised form: 17 May 2012 / Accepted: 21 May 2012 / 
Published: 8 June 2012 



Abstract: This paper is concerned with the reliable finite frequency filter design for 
networked control systems (NCSs) subject to quantization and data missing. Taking into 
account quantization, possible data missing and sensor stuck faults, NCSs are modeled in 
the framework of discrete time-delay switched systems, and the finite frequency I2 gain is 
adopted for the filter design of discrete time-delay switched systems, which is converted into 
a set of linear matrix inequality (LMI) conditions. By the virtues of the derived conditions, 
a procedure of reliable filter synthesis is presented. Further, the filter gains are characterized 
in terms of solutions to a convex optimization problem which can be solved by using the 
semi-definite programme method. Finally, an example is given to illustrate the effectiveness 
of the proposed method. 

Keywords: reliable filtering; networked control system; data missing; quantization; 
sensor faults 



1. Introduction 

In recent years, there has been a growing interest in networked control systems (NCSs), which is a 
class of systems in which sensors, controllers and plants are connected over the network media [1-4]. 
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Due to their advantages such as easy installation, low cost and high utilization, the NCSs have 
widely applications in many application areas, such as manufacturing plants, automobiles and remote 
process, etc. However, these systems require novel control design to account for the presence of network 
in the closed loop, such as network-induced delay (see e.g., [5-8]) and packet loss (see e.g., [9,10]). 
Further, for the NCSs where bandwidth and energy are limited, quantization becomes indispensable. 
Consequently, there has been a lot of researches concerning this issue, (see e.g., [11,12]). 

On the other hand, filtering problem has been playing an important role in control engineering 
and signal processing that has attracted constant research attention, (see e.g., [13-17] and references 
therein). However, it is quite common in practice that measurement outputs of a dynamic system 
contain incomplete observations because of the temporal sensor faults, (see e.g., [18-21] and references 
therein). Therefore, it is natural that the reliable filtering problem in presence of possible sensor faults 
has recently obtained much attention and there have been many results investigating this important issue. 
For example, reliable filtering problems have been thoroughly investigated in [22-24] for linear systems. 
As for nonlinear systems, reliable filtering problems with sensor faults have also attracted many research 
interests [25-27]. 

It should be noted that disturbances considered in those papers are all considered in full frequency 
domain. However, practical industry systems often employ large, complex, or lightweight structures, 
which include finite frequency fundamental vibration modes. Thus, it is more reasonable to design 
reliable filters in finite frequency domain. However, to the best of the authors' knowledge, reliable 
filtering problems for NCSs subject to packet loss and quantization have not been fully investigated, 
especially in finite frequency domain where faults occur frequently. This motivates the investigation of 
this work. 

In response to the above discussions, in this paper, the reliable finite frequency filtering problem for 
NCSs subject to packet loss and quantization is investigated in finite frequency domain against sensor 
stuck faults. Specifically, with consideration of quantization, possible packet losses and possible sensor 
stuck faults, NCSs are modeled in a framework of discrete time-delay switched system. Then, the 
definition of finite frequency I2 gain is given and an analysis condition to capture such a performance for 
discrete time-delay switched system is derived. With the aid of the derived conditions, a reliable filter 
is designed and the conclusions are presented in terms of linear matrix inequalities (LMIs). Finally, an 
example is given to illustrate the effectiveness of the proposed method. 

The reminder of the paper is organized as follows. The problem of system modeling for NCSs with 
packet losses and quantization is presented in Section 2. Section 3 provides sufficient conditions for the 
design of reliable filters. In Section 4, an example is given to illustrate the effectiveness of the proposed 
method. Finally, some conclusions are presented in Section 5. 

Notations: Throughout the paper, the superscript T and —1 stand for, respectively, the transposition 
and the inverse of a matrix; M > 0 means that M is real symmetric and positive definite; / represents the 
identity matrix with compatible dimension; || ■ || denotes the Euclidean norm; P is the probability measure; 
E(-) denotes the expectation operator; I2 denotes the Hilbert space of square integrable functions. In 
block symmetric matrices or long matrix expressions, we use * to represent a term that is induced by 
symmetry; The sum of a square matrix A and its transposition is denoted by He{A) := A + A^. 
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2. System Model and Problem Formulation 

The NCS under consideration is setup in Figure 1 , where the discrete-time plant is of the form: 



x{k + 1) = Ax{k) + Bw{k) 
y{k) = Cx{k) 
z{k) = Ex{k) 



(1) 



where x{k) E R" is the state, y{k) G R™ is the measured output, z{k) E W is the controlled output 
and w{k) E R*^ is the exogenous disturbance which is assumed to belong to /2[0, oo). A, B, C and E 
are known real constant matrices with appropriate dimensions. 

Figure 1. Structure of the Networked Control Systems. 
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In this paper, we make the following assumption: 
Assumption 1. System (1) is stable. 

Remark 2. Assumption 1 is required to get stable dynamics of the filter system. If this assumption is not 
satisfied, a stabilizing output feedback controller is required. 

When the sensors in the NCSs experience faults, we consider the following sensor stuck fault model 
similar to [28], 

y^^{k) = F,y{k) + F,)y,,{k), ^ = 0, 1, 2, . . . , g (2) 
where q is the quantity of the possible fault modes and 

Vsiik) = [ysii{k) ysi2{k) ... ysim{k)f (3) 

with ysij{k){j = 1,2, ... ,m) being the low frequency fault of the kth sensor. Further, is defined as 

Fi = diag {Fn, Fi2, . . . , Fj^} 
Fik = 0 or 1, k = 1,2, . . . ,m 

It is also assumed that, as shown in Figure 1, the measurement signals will be quantized before 
transmitting via the networks wherein data missing may occur. The following logarithmic quantizer as 
proposed in [29] is applied, 

Uvo if T^/^o < V < Y^/vo 
q{v) = < 0 if = 0 (5) 

I -q{v) if i; < 0 
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where the parameter 0 < p < 1 is termed as quantization density and 

5, = (l-p)/(l + p) 

From [29], we can obtain 

q{v) 



(6) 



(7) 



where Ag G [— 5q, is a suitable model for the logarithmic quantizer q{v) with parameter 5q. 

Therefore, the faulty measurements together with quantization and the data transmission in the 
networks can be described by 

(A;) = a{I + A,)/'(A;) + (!-«)(/ + A,)/'(A; - 1) (8) 

where a G R is a Bernoulh distributed white sequence with 

P(a = 1) = E(a) = 5 

P(a = 0) = 1 - E(a) = 1-5 

Specifically, if a = 1, the quantized signal (/ + Aq)y^^{k) is successfully transmitted, otherwise the 
transmission fails, i.e., the phenomenon of data missing. 

Remark 3. The description of data transmission (8) was introduced in [30]. It can be seen that the 
output y{k) of the system model is (/ + Aq)y^* (k) with probability 6 at k-th sampling time, and the value 
(/ + Aq)y^^ [k — 1) with probability 1 — 5. Obviously, if the binary stochastic variable a takes the value 
0 consecutively at different sample times, the consecutive data missing would occur 

In this paper, the following reliable filter is constructed: 

x{k + l) = Afx{k) + Bfy^^{k) 



i{k) = Cfxik) 



(9) 



where Af, Bf and Cj are filter parameters to be designed. 

Denoting ({k) = [x'^{k) and e{k) = z{k) — z{k), then the filtering error system for the 

zth fault mode can be described by the following two subsystems. 

5*1: No packet dropout occurs. 

C{k + 1) = AuC{k) + Au^C{k - 1) + B^w{k) + Bsiys^{k) 

e{k) = cm 

S2'. Packet dropout occurs. 

C{k + 1) = A2iC{k) + A2diC{k - 1) + B^w{k) + Bs^ys^{k) 
eik)=CC{k) 



where 



An Aidi A2i A 



2di 



A 0 
Bf{I + Aq)F,C Af 



0 


0 


A 


0 


0 


0 


0 
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Due to packet drop-out, the filtering error system can be seen as combined by subsystem Si and 5*2, 
wliich can be lumped into the following discrete time-delay switched system: 



C(/c + 1) = A^^,C{k) + A^^d^C{k - 1) + B^w{k) + BsiVsiik) 
e{k) = CC{k) 



(10) 



where af, is switching signal with ak G S = {1,2} being a piecewise constant function. 

Next, we will discuss how to design the filter parameters Af, Bf and Cf. In order to formulate the 
problem clearly, the following definitions are first given. 

Definition 4. (Asymptotical Stable) System (10) is said to be asymptotical stable under switching signal 
ak, if the solution satisfies 

lim||C(fc)|| =0 (11) 

Definition 5. Ding et al. [17] Let 'j > 0 be a given constant, then the filtering error system (10) is said 
to have a finite-frequency I2 gain 7, if inequality 



(12) 



A:=0 



fe=0 



holds for all solutions of Equation (10) with w{k) G I2 such that the following hold 
(i) For the low-frequency range \6\ < -di 

k=0 



j^icik + 1) - amcik + 1) - c(^))^ < {2s^n^r Yl 

( ii) For the middle-frequency range < 0 < ■§2 

00 

eo^^ J2(C(k + 1) - e^''^C(^))(C(A: + 1) - e-^^^ak)f < 0 



(13) 



(14) 



k=0 



where 1^^ = {^2 - ^i)/^. 
( Hi) For the high-frequency range \6\ >'dh 



(15) 



Y,m + 1) - c(^))(c(^ + 1) - mf > {2s^nfr y1 a^n^f 

k=0 k=0 

Now, the reliable filtering problem to be addressed in this paper can be formulated as follows: 
Design a stable reliable filter (9) such that, for the quantization error, possible data missing and sensor 

faults, the filtering error system (10) is asymptotical stable, and with a prescribed finite-frequency I2 gain 

7i from w(^k) to e(^k) by satisfying the following specification 



E 11^(^)112 <7iE 11^(^)112' ^^'^ 



(16) 



k=0 



k=0 
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and a prescribed low frequency I2 gain 72 from Usiik) to e{k) by satisfying 

00 00 

$^||e(A;)||2<72'Ell^-(^)ll2' V |^^.| < ^ 



(17) 



fc=0 



fc=0 



where 6^ and 6^^ represent the frequency of the disturbance and the stuck faults, respectively. 

Before ending this section, the following lemmas will be first given to help us to prove our 
main results. 

Lemma 6. (Finsler's Lemma) For x G R",£ G R"^^",!! G R"^™, let il^ be any matrix such that 
il-'-il = 0. The following statements are equivalent: 

(i) < 0, Vil^x = 0, X ^ a 

(ii) il^m^^ < 0, 

(Hi) 3/i G R" : £ - /iilil^ < 0, 

(iv) 32) G R'"^" : £ + il2) + 2)^il^ < 0. 

Lemma 7. Given the matrices E and F with appropriate dimensions, then 

EAF + F^A^E^ < 0 
where AA^ < /, if and only if there exist a scalar e > 0 such that 



e—2F 



< 0 



3. Main Results 

In this section, the reliable filtering problem proposed in the above section will be investigated. 

Lemma 8. Consider system ( 10) for i = 0,1, . . . ,q and a given scalar 71 < 0, then Equation (16) holds, 
i.e., system (10) is with a finite frequency I2 gain 71, if there exist matrices Va^i = 'Paki' ^« ~ ^'^'^ 
Qaki = Qau > 0, cTfc G {1, 2} such that the following inequalities hold 



10 0 



10 0 



+ 





" c 


0 


0 


+ 










0 


0 


I 



n 



" / 


0 


0 


T 




0 


0 


/ 


0 




0 





c 


0 


0 


0 


0 


I 


I 


0 


0 


0 


I 


0 



(18) 



< 0 



where 11 



I 

0 



0 

-ill 



and 



(i) For the low-frequency range |6'^| < 'dyji 

'V, 



Qo 



(19) 
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( ii) For the middle-frequency range < 9^ < ^2 



-V, 



Q 



(20) 



where A = {A + ^1)/% K = (??2 - ^?i)/2. 
(in) For the high-frequency range > '&wh 



-V. 



o 



(21) 



Proof. We first consider tlie middle-frequency case for system (10) with ysi{k) = 0. Assume 



Equation (14) holds, pre- and post-multiplying it by 



and its transpose, we can derive 



e{k)iv,,, + R,)ak) - e{k + i)p.,^,.c(A: + 1) - e{k - i)Riak - 1) 

+ e^{k)e{k) - -ilw'^{k)w{k) 

+ triQ^^^ie^'^^mCik + 1) - e~^'>X{k + 1)C^(A;) - 2cosd^ak)C {k))} < 0 

Summing up Equation (22) from 0 to 00 with respect to k, it is easy to obtain 

00 

{^{e^{k)e{k) - ^lw^{k)w{k))) + tr{Q,^iZ} < 0 

A:=0 

since system (10) is asymptotical stable and ^(0) = 0, where 

00 

Z := Y.{e^'X{k)e {k + 1) + e-^^^^^k + 1)C^(A;) - 2cos^^a^)Cm 

k=0 



(22) 



(23) 



(24) 



It is easy to prove that —Z is equal to the left-hand side of Equation (14), thus we have Z > 0. Further, 
from Qa-f^i > 0, one can obtain that the term tr{Q„i,iZ} > 0 while Equation (14) is satisfied. Hence we 
have YlkLo (E(e^(A;)e(A;) — '^\w'^ {k)w{k)) < 0, which is equivalent to the condition Equation (16) for 
middle-frequency in Definition 5. 

Similarly, by choosing "di := —i^wi, ^2 '■= "dwi for low-frequency case and di := 'dwh, ^2 '■= 2tc — d^h 
for high-frequency case, respectively, the results for these two cases can be derived by following the 
same procedure of the above proof. This completes the proof. □ 

Remark 9. Lemma 8 presents an analysis condition for finite frequency I2 gain of system ( 10), where less 
conservatism is introduced compared with the existing full frequency conditions when frequency ranges 
of disturbances are known. 

Remark 10. The sufficient condition, which guarantees a prescribed low frequency I2 gain from ysi{k) 
to e{k) for system (10), can be obtained by following the same process of Lemma 8 and utilizing relevant 
system matrices. 
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3.1. Finite Frequency Performance from w{k) to e{k) 



In this section, sufficient conditions to capture the finite frequency performance from w{k) to e{k) for 
system (10) will be derived with the aid of Lemma 8. 

Theorem 11. Consider system (10) in fault free and faulty cases (i.e., i = 0, 1, . . . ,q) for given low 
frequency range \9\ < 'dwi, which is with a prescribed I2 gain •yifrom w{k) to e{k), i.e., the condition 
(16) holds if there exist a scalar £1 > 0, matrices X, y, Af, Af, Bf, Cf and 



V = V 



ni = Ui 





* 








* 







V = V 





* 


^Ki2 




QkH 


* 







> 0 



with L, K E {1, 2} such that the following conditions hold 

^ < 0 

where 
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* 


* 
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QkH — ^ 
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Qkvi - y 
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^ = 
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* 




0 
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-7^^2 


-7^^3 
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* 
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0 




BTyT 
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0 


-ill 
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* 




0 


0 


E 


-Cf 
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0 


0 


-/ 


* 




0 


0 






0 


0 


0 


0 


-eil 



and 





K 


= 1) 


= V^a - 2cosdu,iQ.ii + Ua + S^eiC^F^^FiC + He{XA + BfF^C), 


Vl/3( 


K, 


= 2) 


= V^ii - 2cosdu,iQ.ii + 7^a + He{XA), 




K, 


= 1) 


= V^i2 - 2cosd^iQ^i2 + 7^.2 + BfFiC + ^ J + yA, 




K, 


= 2) 


= V^i2 - 2cosi^^iQ^i2 + 7^,2 + ^/ + yA 








- 2cosi^iQ^i3 + 7^i3 + He{Af), 


m 


K 


= 1) 


= nlin = 1) = 0, vl;3(^ = 2) = = 2)= C^F^Bj 


m 


K 


= 1) 


= -TZa, = 2) = -7^,,l + 5]eiC^F^FiC 



(25) 



Proof. It is shown, from Lemma 8, that the condition (16) can be reached if Equation (18) holds. Further, 
the inequality (18) can be rewritten to 



j^r < 0 



(26) 
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where J 





I 


0 


0 




0 


/ 


0 




0 


0 


/ 



A 



and 



Ski 

0 
0 



* * 
7r/ 



2cos^^iQ^, + ni + C'C * 
0 0 



-ill 



Exploiting Lemma 6 and explicit null space bases calculations on Equation (26), it is easy to get that 
Equation (26) holds if and only if 



A + He(^J^[0 0 0]) < 0 



(27) 



that is 



Q.^ - W V^^- 2cosd^iQ,i + 7^, + C^C + He{WA^^ 
0 



* 

0 



-ill 



< 0 



where W is a matrix variable introduced by Lemma 6 and J-*- = 
Performing Schur's complement on Equation (27) yields to that 

-Va * 
0 BlW^ 

0 c 



I A^.i B-,,, A, 



ndi 



is utilized. 



0 
0 



* 

0 



< 0 



Partitioning W as the following form 



y M 



(28) 



defining the following new variables 



Af=MA 



Bf=MB 



f 



and applying Lemma 7 on Equations (27), (25) can be reached easily. This proof is completed. 



(29) 

□ 



Remark 12. In Theorem 11, by introducing a variable W, the coupling between the Lyapunov matrices 
and the filter gains will be eliminated, which does not present any structural constraint. 

The previous Theorem 11 presented the condition to capture the low frequency performance. 
Similarly, conditions for middle frequency and high frequency performance are presented in the 
following two corollaries. 
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Corollary 13. Consider system (10) in fault free and faulty cases (i.e., i = 0,1, . . . , q)for given middle 
frequency range < \9w\ < {^2, which is with a prescribed I2 gain •yi from w{k) to e{k), i.e., the 
condition (16) holds if there exist a scalar e\ > 0, matrices X, y, M, Af, Bf, Cj and 



V = V 
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7^^2 





V, 
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QkA — Qui 
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QkiS 



> 0 



with L, K G {1,2} such that the following conditions hold 

^ < 0 



where 



and dr 
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-eil 



1) = ^Kil - ^COS^uiQtiil + Tin 

2) = V^ii - 2cosi^^Q^ii + n^l + He{XA), 

1) = V^i2 - 2cosdM.i2 + ni2 + BfFiC + A] 

2) = V^i2 - 2cos^M^^2 + 7^,2 + ^/ + yA, 



5leiC^F^F,C + He{XA + BfF.C), 



V|/4 



C^F^'BJ 



^lin = 2) = nlin = 2) 
1) = -7^.l, = 2) = -7^a + Sle^C' Ff F^C 

Proof. By following the same lines of Theorem 1, it is immediate. 



(30) 



□ 



Corollary 14. Consider system (10) in fault-free and faulty cases (i.e., i = 0, 1, ... , q) for given high 
frequency range |6'^| > which is with a prescribed I2 gain •yifrom w{k) to e{k), i.e., the condition 
(16) holds if there exist a scalar ei > 0, matrices X, y, Af, Af, Bf, Cf and 
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with L, K E {1, 2} such that the following conditions hold 



^ < 0 



(31) 
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where 



and 
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1) 

2) 

1) = P«2 + 2cos^^hQ,i2 + 7^i2 + BfFiC + ^ J + yA, 

A] + yA, 



V^a + 2cosKhQnii + 7^il + SleiC^Fj^FiC + He{XA + S/F^C), 
V^ii + 2cosd^hQnii + 7^^l + He{XA), 

2) =V^i2 + 2cOS^^hQKi2+ni2 

V,i3 - 2cosm.i3 + T^i3 + He{Af), 
<fl{K = 1) = nl{K = 1) = 0, mn = 2) = nlin = 2) 

^^(/s: = 1) = -7^a, = 2) = -7^,;l + dleiC'F^F.C 

Proof. By following the same lines of Theorem 1, it is immediate. 



C^F^BJ 



□ 



3.2. Low Frequency Performance from ysi{k) to e{k) 



In this subsection, sufficient conditions to capture the low frequency performance (17) for system (10) 
will be deduced. 

Theorem 15. Consider system (10) in faulty cases (i.e., i = 1, . . . ,q) for given low frequency range 
\0s\ < ^si, which is with a prescribed low frequency I2 gain 72 /or nonzero ysi{k), i.e., condition (17) 
holds, if there exist a scalar 62 > 0, matrices X, y, Af, Af, Bf, Cf and 
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(32) 
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and 



fill 



1) = 

1) = 

1) = 
K = 2) = 

1) = 



VsKii - 2cos^siQsKii + Tlsii + He{XA), 

Vsni2 - 2cos§,iQsKi2 + 7^.i2 + BfFiC + A) + yA, 

0, fi^c^^2l 



^1 = 



'- 1) — —T^sil, 

-7i/ + &(/-F,f(/-F.), 



-7^ 



sil + Sq£2C Fj FjC, 



Proof. It is easily derived, from Lemma 8 and Remark 10, that the condition (17) holds if 



A. 



■CTfcJ 
/ 



A(jf,(ii 
0 



A(jf^i A(jf,iii Bgi 


+ 


c 


0 


0 


± 

n 


c 


0 


0 


10 0 




0 


0 


I 


0 


0 


I 




1 T 

















I 

0 



0 

I 



0 



0 



I 

0 



(33) 



< 0 



where 'E^ has the similar structure as in Lemma 8 for low frequency case. 

Following the same process in Theorem 1 1, we first rewrite the inequality (33) to the following form. 



£e£^ < 0 



(34) 



where £ 



Af^i 



e 



/ 0 0 
0/0 
0 0/ 



and 



-V 

' sn 
QsKi 

0 
0 



QsKi * * 

2cos^siQsKi + nsi + C^C * * 

0 -TZsi * 

0 0 -7|/ 



Exploiting Lemma 6 and explicit null space bases calculations on it, we have Equation (34) is 
equivalent to 

S + //e(^£^[0 0 0])<0 (35) 

that is 



-Vsd * 

0 
0 
0 



AL;w^ 



* 

0 
0 



-ill 

0 



* 



< 0 



/ Ani A, 



ndi 



B 



fi 



is utilized. 



where W is a matrix variable introduced by Lemma 8 and £ = 

By applying Lemma 7 on Equation (35), the inequality Equation (32) can be reached easily. This 
proof is completed. □ 
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3.3. Stability Condition 



Since Theorems 1 and 2 presented in the above section cannot guarantee the stability of the system 
(10), in this subsection, asymptotical stability conditions for system (10) will be presented. 

Theorem 16. Consider system (10) in fault free and faulty cases (i.e., i = 0,1, ... ,q), it is asymptotical 
stable when w{k) = 0 andysiik) = 0, if there exist a scalar > 0, matrices X, y, Af, Aj, Bf, Cf and 



V. 



T 
an 



V. 



T 



TZai 



'Paul 


* 


'Pad2 






* 








* 


TZai2 





> 0 



> 0 



> 0 



with L, K E {1, 2} such that the following conditions hold 

$ < 0 

where 



(36) 



V, 



aiil 



He{X) 



A- 



* 


* 


* 


* 


* 


* 


- He{M) 


* 


* 


* 


* 


* 






* 


* 


* 


* 




—VaK%2 + T^aO. 




* 


* 


* 




0 


0 




* 


* 


0 


0 


0 






* 




0 


0 


0 


0 


-83! 



and 





[k, = 


1) 


= A^A"^ + C^F^BJ, 




= 2) 


= A^X^ 






[k, = 


1) 


= A^'y^ + C^FfBj, 




= 2) 


= A^y^, 






[k, = 


1) 


= -Va.^l + 7^„l + dle^C^F^F.C, 




= 2) 








[k = 


1) 






= 2) 


= ^li^ = 


-- 2) = C^i^^Sj, 




[k = 


1) 


Traill 


^1{k 


= 2) 


= —T^ail 


+ 5le,C^FrF,C. 



Proof. Consider the following Lyapunov functional candidate for w{k;) = 0 and ysiik) = 0 

V,,{k) = e{k)Va.^m + ^{k - l)na^ak - 1) 



(37) 



The forward difference of the Lyapunov functional Vi{k) along the trajectories of the system (10) is 
given by 



AV^,{k) = V,,{k + l)-V^iik) 

= e{k){Al,VanA,^ - Van + 7^..)C(A:) + {k){Al,VanA^di)ak " 1) 

+ e{k - l)iAl,,VanA,i)ak) + e{k - l){Al,,VanA,,, - 7^,,)C(A; - 1) 

= e{k)rMk) 



(38) 
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where e(A;)= (^{k) (^{k-l) 



and 



T 



, V A ^ —TZ 



By following an opposite direction to the proof for Theorem 1 1 and exploiting Schur's complement 
and Lemma 7 on Equation (36), we have T^j < 0, which implies that AVKi(fc) < 0. Therefore, from 
Equation (38), one can easily obtain, for a sufficiently small scalar p > 0 and C,{k) ^ 0, that 



AK.(fc)<-p||C(A;)|p 



(39) 



Then, from Definition 4, the asymptotical stability of system (10) can be established. This proof is 
completed. □ 

3.4. Algorithm 

Based on the above analysis, a set of optimal solutions Af, Bf and Cj can be obtained by solving the 
following optimization problem for given 5g: 



min a'-fi + 672 

s.t. (25), for 2 = 0, 1, . . . , g, 

(32), fori = 1,2,. .. ,g, 

(36), fori = 0,1,..., g. 



(40) 



where a and b are weighting factors. 

Then the dynamic output feedback controller gains can be computed by the following equalities: 



Af = Af-^Af, Bf = Af-^Bf, C 



f 



C 



f 



(41) 



Remark 17. On the other hand, we can obtain a coarser quantizer through solving the following 
optimization problem for given finite-frequency I2 gains 71 and 72 



max 



s.t. (25), /ori = 0,1,..., g 
{32), for i = 1,2,..., q 
{36), fori = {),!,..., q 



(42) 



4. Example 



In this section, an application and simulations are given to illustrate the effectiveness of the 
proposed methods. 
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The utilized model is the F-404 aircraft engine described by the following state space model [31], 





-1.4600 


0 


2.4280 




-0.09 


x{t) = 


0.1643 


-0.4000 


-0.3788 


x(t) + 


-0.14 




0.3107 


0 


-2.2300 




0.02 


y{t) = 


1 1 0 
0 1 1 


x{t) 









w{t) 



(43) 



z{t) 



0.2 -0.7 
0.3 0.6 



0.5 
-0.4 



x{t) 



Assume the sampling period is /i = Is, and packet transmission is as shown in Figure 2, which is 
subject to the rate of packet lost 1 — 5 = 0.32. 

Figure 2. The number of the lost packet. 




20 25 30 

time steps (k) 



50 



For given = 0.2, = 0 and 5q = 0.6, solving the optimization problem (40), we can obtain the 
optimal value for low frequency performances are, respectively, 71 = 0.0825 and 72 = 0.0378 with the 
corresponding reliable filter parameters 



x{k + l] 



z{k) 



0.1454 -0.0418 0.0447 
0.7084 0.1740 0.1958 
0.0246 0.0061 -0.0620 

0.7055 0.8941 -0.0583 
-1.2318 -0.7672 1.4227 



x{k) + 
x{k) 



0.0813 0.0813 
0.2741 0.2741 
-0.0049 -0.0049 



(44) 



In the following, the system is simulated in low frequency domain, where the faults always occur, 
under the following two fault modes with zero initial condition and the disturbance input w{k) is 



w{k) 



0.5sin{k), 
0, 



10 < A; < 15 
otherwise 



which is shown in Figure 3. 
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Figure 3. The disturbance input w{k). 




Mode 1. The first sensor being stuck at 0, i.e., 

F = 

and ysi{k) = 0 when k > 30. 
Mode 2. The second sensor being stuck at 0, i.e.. 



0 0 
0 1 



1 0 
0 0 



and ysi{k) = 0 when k > 30. 



10 20 30 40 

k 



50 



The controlled outputs and the corresponding estimations for both the two fault modes are shown in 
Figures 4 and 5, respectively, where the blue solid lines are the controlled output while the red dashed 
lines are their estimations. 

It is easily seen from these figures that all the expected system performance requirement are well 
achieved, which shows the effectiveness of the proposed method. 
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Figure 4. The z{k) (blue solid) and z{k) (red dashed) when faulty mode 1 occurred. 
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Figure 5. The z{k) (blue solid) and z{k) (red dashed) when faulty mode 2 occurred. 
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5. Conclusions 

In this paper, the reliable finite frequency filtering problem for NCSs subject to quantization and 
packet have been studied with finite frequency specifications when sensor fault would occur. The 
considered NCSs have been first modeled as a discrete time-delay switched system. Subsequently, a 
sufficient condition to characterize the finite frequency I2 gain has been presented. Then by virtues of 
the derived condition, a procedure of reliable filter synthesis has been derived in terms of LMIs. Finally, 
an example has been given to illustrated the effectiveness of the proposed method. 
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